"Bad Metal" Conductivity of Hard Core Bosons 
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Two dimensional hard core bosons suffer strong scattering in the high temperature resistive state 
at half filling. The dynamical conductivity u(cj) is calculated using non perturbative tools such 
as continued fractions, series expansions and exact diagonalization. We find a large temperature 
range with linearly increasing resistivity and broad dynamical conductivity, signaling a breakdown 
of Boltzmann-Drude quasiparticle transport theory. At zero temperature, a high frequency peak in 
ct(oj) appears above a "Higgs mass" gap, and corresponds to order parameter magnitude fluctuations. 
We discuss the apparent similarity between conductivity of hard core bosons and phenomenological 
characteristics of cuprates, including the universal scaling of Homes et al. (Nature 430, 539 (2004)). 

PACS numbers: 05.30. Jp, 72.10.Bg,74.72.-h 



I. INTRODUCTION 



Quantum transport in condensed matter is largely based 
on the paradigms of Fermi and Bose gases. Boltzmann 
equation for the conductivity is valid in the weak scat- 
tering regime where it yields a Drude form [l|, Q , 



m* 1 — hot 



(1) 



where T is temperature, uj is frequency, and q, n, m* , t(T) 
are charge, density, effective mass and scattering time of 
the constituent quasiparticles. 

Interacting bosons in a strong periodic potential 
may suffer strong enough scattering which invalidates 
Boltzmann-Drude theory. An example is provided by 
the two dimensional Hard Core Bosons (HCB) model at 
half filling. While it is established that the ground state 
is a bone fide superconductor [H, 0, [H, 0, 0| , the resistive 
(normal) phase involves strongly interacting bosons and 
vortex pairs. Previous work [H, Q showed that the lat- 
tice (umklapp) scattering dramatically increases vortex 
mobility. At half filling, it produces an abrupt reversal 
of the Hall conductivity, and doublet degeneracies asso- 
ciated with each vortex. 

HCB models may be experimentally relevant to cold 
atoms in optical lattices (ToL Hl| . underdoped cuprate 
superconductor s fl2. [T3 , Il4|. low capacitance Josephson 
junction arrays 15l . 16l|. and disordered superconducting 
films An added advantage of HCB, is that they are 
described by a quantum spin-half XY model which is 
amenable to tools of quantum magnetism. 
It is the purpose of this paper to compute the conduc- 
tivity of HCB at half filling. We apply and test a set of 
non perturbative approaches, including continued frac- 
tion representation, series expansions and exact diago- 
nalization. By studying the dynamical structure of the 
Kubo formula, we can construct well-converging approxi- 
mants which agree with high accuracy sum rules in a wide 
regime of temperature. The conductivity at high temper- 
atures is obtained to order of T~ 3 . Near the supercon- 
ducting transition, it is matched with the critical conduc- 
tivity which was derived by Halperin and Nelson [l8[ . At 
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FIG. 1: Temperature dependent resistivity of two dimensional 
Hard Core Bosons at half filling. High temperature line (blue 
online) is calculated up to order 1/T, with error margins de- 
picted by dashed lines. The critical region above the BKT 
transition, (red online) is given by the vortex plasma the- 
ory of Halperin and Nelson (HN). For illustration, a layered 
system with weak interlayer coupling shows a rapid rise in re- 
sistivity above the three dimensional transition temperature 
T^ D , as depicted by a solid black line. 



zero temperature, the dynamical conductivity is obtained 
from the relativistic Gross-Pitacvskii field theory, and a 
variational fit to 12th order moments. 
Our key results pertain to the qualitative effects of strong 
scattering on the conductivity. At high temperatures, the 
conductivity goes as 



cr(T» w 0.91 



q 2 tanh(fiu/(2T)) 



■ cxp 



(2) 



where f2 is a high frequency scale. Eq. when fit 

at low frequencies and temperatures, to the form (fTJ), 
yields an "effective scattering rate" which is equal to 2T. 
The resistivity as shown in Fig. [T] exhibits "bad metal" 
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fl9| behavior: it exceeds the boson Ioffe-Regel limit [2(| 
of R > h/q 2 , without saturation [2l], [52]. The resistiv- 
ity's linear slope defines the proportionality coefficient 
between the zero temperature superfluid stiffness, and 
product of T c and the normal state conductivity near T c . 
We note that linearly increasing resistivity [231 ] and an 
analogous scaling of superfluid stiffness with conductiv- 
ity called "Homes law" 24|, have been observed in the 



cuprate family, as depicted in Fig. [9l 
In addition, at zero temperature we find a small conduc- 
tivity peak above a "Higgs mass" gap. This peak is as- 
sociated with order parameter magnitiude fluctuations. 
These are analogous to coherence oscillations observed 
in cold atoms during rapid Mott to superfluid quenches 
[251 . 12H I27I ] . We speculate that perhaps the mid infra-red 
peak observed in some cuprates at low temperatures [28| , 
might arise from these magnitude fluctuations. 
However we emphasize that similarities between HCB 
and cuprates are only suggestive. This work does not 
include effects of fermion quasiparticles, magnetic exci- 
tations, and inhomogencities which are experimentally 
important. 

This paper is organized as follows: Section ITT1 introduces 
the HCB model and lists some of the established ther- 
modynamics. Section IIIII discusses the Kubo formula in 
general, and derives the moment expansion, continued 
fraction recurrents, and orthogonal polynomials which 
can be used to evaluate a non perturbatively and gen- 
erate a high temperature expansion. Section IIVI derives 
the particular recurrents for the HCB at half filling. Jus- 
tification for the rapidly converging harmonic oscillator 
expansion is provided. Sections IIIIIIVI are quite techni- 
cal, and could be avoided at first reading. The results 
of our calculations are provided in the following sections. 
Section |V] plots the dynamical conductivity at high tem- 
peratures, and explains three approaches which converge 
to the same curve. Section |VT1 obtains the resistivity as a 
function of temperature. Section TVIII obtains the dynam- 
ical conductivity and "Higgs mass" at zero temperature. 
Section IVIIII summarizes the key results and discusses 
their possible relevance to the dynamical and DC con- 
ductivity of cuprate superconductors. 



II. HARD CORE BOSONS 



A minimal model of HCB hopping with Josephson cou- 
pling J, coupled to an electromagnetic field Aij, is the 
gauged S = \ quantum XY model, 



H = -2J£ 



a\aj + H.c. 



(i,0) 



2.1 (<" iA ' ! I-''-) 

(id) 



(5) 



where q is the boson charge (= 2e, for electronic su- 
perconductors) and we use units of h = c = I. Here we 
consider (i, j) to be nearest neighbor bonds on the square 
lattice. In the absence of a chemical potential (i.e. a Zee- 
man field coupled to ^ 5?), the Hamiltonian describes 
a density of half filling (zero magnetization) , with half a 
boson per site. 

The uniform current operator, for A = 0, is given by 



Jx 



2iqJ 



AqJ 



E 

(id) 

E (#5 



ajdj — a^a, 



y _ qy qx \ 
r+x °r°r+x) 



(G) 



We note that in one dimension, the current of HCB is a 
conserved operator since, 



[H W ,J X ] =0, 



(7) 



and hence the real conductivity trivially vanishes at all 
finite frequencies and temperatures. In two and higher 
dimensions, this is not the case: the conductivity has non 
trivial dynamical structure. 

Below we review some established results for the thermo- 
dynamic properties of the two dimensional quantum XY 
model which are relevant to the conductivity. 



A. Superfluid stiffness 

It is widely believed that at zero temperature H has long 
range order. Thus, at low temperatures T > 0, the two 
dimensional boson superfluid stiffness p s , which has units 
of energy, is finite: 



Hard core bosons are defined on lattice sites i = I, ... TV 
with restricted occupation numbers, m = 0, 1. The con- 
strained creation operators a\ are represented by spin 
half raising operator 



at 



s; 



Thus, their commutation relations are 

[S i ,Sj\ = —5ij2S* 
hi, a] = 6ij (1 - 2n.j) 



(3) 



(4) 



p s = q 



(dA x f 



>0, 



(8) 



A=0 



where F is the free energy in the presence of a uniform 
field A = A x ~k. A classical (mean field) approximation 
at zero temperature yields a non-monotonous density de- 
pendence, 



pf(0,n) =4Jn(l-n). 



(9) 



where n is the mean boson occupation (filling). Half 
filling n = i is ''optimal", with maximal p s . Quantum 
corrections to pf(0, ^) enhance it by about 7% [1, 0]- 
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B. BKT Transition 

The static order parameter correlations of ([5]) are de- 
scribed by a renormalized classical XY model. At low 
temperatures, correlations decay as a power-law in dis- 
tance. Ps(T) decreases with T, until it falls discon- 
tinuously to zero at Tbkt, the Berezinskii-Kosterlitz- 
Thouless (BKT) transition temperature 0, [3(| [3l[ . At 
half filling, quantum Monte-Carlo (QMC) simulations of 
(J5]) have determined 0, [|[ , 

T BKT ~ 1.41 J. (10) 

Just below the transition, a universal relation holds: 



Ps(Tbkt) = —Tbkt- 



(11) 



field theory pjj, and the mean field superfluid stiffness 
goes as 



Ps 



AJn 



(17) 



where mj, a are the boson effective mass and lattice con- 
stant respectively. Similarly, by particle hole transforma- 
tion l|15p , H simplifies into a model of weakly interacting 
holes asn-^ 1. 

ft is around half filling, however, that higher order inter- 
actions in p6[) and lattice effects become relevant. The 
Galilean invariant (non-relativistic) GP theory fails to ac- 
count for important umklapp scattering processes which 
lead to Hall effect cancellation, and doublet degeneracies 
of vortex states H . 



For T > Tbkt, Ps = 0, and the correlation length is 



£ ~ A exp 



B 



(12) 



y/(T — Tbkt) /Tbkt J 

where A = 0.285 and B = 1.92 as determined by QMC 
1- 

In multilayered systems with weak intcrlaycr cou- 
pling [13, HU, [33| J c « J, the 3D transition temperature 
is higher than Tbkt by a factor, 



T 



3D 



Tbkt 1 + 



B 2 



log 2 (0.144 J/J c ) 



(13) 



C. Boson particle- hole symmetry 

The charge conjugation operator 

transforms particles into oppositely charged holes: 



n% -* (1-rii), 
H[qA,n] -► H[—qA, 1 -n}. 



(15) 



It follows therefore that under reflection about half filling, 
the Hall conductivity reverses sign, while the superfluid 
stiffness and longitudinal conductivity are invariant. 
In the low density limit n « 4, the HCB are effectively 
unconstrained [34],[35| as seen by (J2J). This can be demon- 
strated by expanding the Holstein-Primakoff bosons rep- 
resentation of spin half, 



1-n? 



b\(l 



0(n 



b\2\ 



(16) 



Truncating the expansion (|16p . and inserting it into H, 
turns it into an interacting soft bosons model. At low 
densities, the low excitations of H do not feel the lat- 
tice. Thus, the long wavelength properties are well de- 
scribed by the Galilican invariant Gross-Pitacvskii (GP) 



III. NON PERTURB ATIVE ANALYSIS OF THE 
KUBO FORMULA 

A. Current fluctuations function 

Thermodynamic averages are denoted by 

((•)>„ = Itt [e-^(-)] , Z = Tre-(> H , (18) 

where (3 = l/T. The linear response Kubo formula for 
the longitudinal dynamical conductivity is [3(| 



(-q 2 K x )p - A xx (P,u; + ie) 



u> + ie 



where 



4 7*- 



N 



v \ 

r+x/ ' 



(19) 



(20) 



is the x-kinetic energy, and A is the retarded current- 
current response function 



A xx (z) = / dte izt ([J x (t),J x (O)]) 

JQ 

J x (t) = e iHt J x e- iHt , (21) 

where the HCB current operator is given by 
The (real) conductivity is given by 

a(M = q \pM5(") + t&nHM2) G"(f3,u;) 

UJ 

poo 

G"{p,oj) = \ \ dte-^({J x (t),J x (0)})p. (22) 



where {-, ■} is an anticommutator. 

The current fluctuations function G"(j3,uj) of HCB at 
half filling, is our primary object of attention. We cannot 
rely on a nearly free quasiparticle basis about which to 
expand Eqs. (|5|) or ((6|). We therefore turn our attention 
to non perturbative approaches. 
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B. Moments expansion 

It is advantageous in our case, to analyse the dynamical 
structure of G"(u>) in the Operator Hilbcrt Space (OHS). 
The OHS is a linear space of HCB (spin half) operators, 
denoted by capital roman letters, A, B, . . ., which are the 
" hyperstates" of the OHS. In this paper, we use two dif- 
ferent inner products, 
(i) The infinite temperature product, 



1 



(AB) 00 = w Tt(A^B) 

(ii) The zero temperature product. 

(A,B) o = (0\A^B\0), 



(23) 



(24) 



where |0) is the ground state of the Hamiltonian H . It 
is easy to verify that both definitions obey the Hilbcrt 
space conditions for an inner product. Henceforth, wc 
unify the notations, and drop the subscripts 0, oo. 
''Hyperoperators", denoted by capital script letters, are 
linear operators which act on hyperstates of the OHS. 
The Liouvillian hyperoperator C is defined by its action 
on any hyperstate A as, 



LA = [H, A] 



(25) 



By hermiticity of H, and the cyclic property of the trace, 
the Liouvillian C is Hermitian for both definitions of the 
inner products (|23l24p . 



{A, LB) = (B,CA)* 



(26) 



Therefore, £ has a real eigenspectrum. 
The time dependent current operator in the Heiscnbcrg 
representation J x (t), is compactly expressed using the 
evolution hyperoperator, 



Jx(t) 

The hyper- resolvent Q(z), 



e lCt J, 



1 



is related to the evolution hyper-operator by 
ACt f dz 



(27) 



(28) 



(29) 



The contour surrounds the spectrum of C, which by Eq. 
(126p lies on the real axis. 

The complexified current fluctuations function is 

G(/3, z) = ilr {e-> 3H {J x ,G(z)J x }) . (30) 

Eq. (|22|) is recovered by its imaginary part on the real 
axis, 



G"{p,w) 



-ilmTrfe-^jj,, i— 

Z \ w — L + le 



Jx 



(31) 



A direct \/ui expansion of (u> — £) _1 does not yield the 
required imaginary function [37j. To extract G" one uses 
complex analysis: 

Idzz k Q(z) = I— z k jr (-Y =2niC k , (32) 
^ j z n=0 \ z / 

and take the contour around the real axis, to obtain the 
sum rules for all k = 0, 2, 4 . . . oo, 



00 dw 



uj k G"((3,Lu) = ({j x ,C k J x })p.=pi k {{3). (33) 



All odd-fc moments vanish by symmetry of G"(uj). The 
moments are static (equal time) correlators, which 
can be evaluated numerically or by series expansions. It 
is possible, in general, to compute only a finite number of 
/Xfe's. The remaining task is to derive converging approx- 
imants to G"{(3,uj) by, in a sense, inverting Eq. 



The Liouvillian Matrix 



In order to invert (|33|) . we use the structure of the Li- 
ouvillian matrix. A tridiagonal matrix representation of 
the Liouvillian is constructed as follows. We define the 
root hyperstate by the current operator, 



O = 



1 



\J (Jx, Jx) 



Jx 



(34) 



An orthonormal set of hyperstates O n can be generated 
by sequentially applying C and orthonormalizing by the 
Gram-Schmidt procedure, 



O n+ i = c n+ i [CO n - A„0„_i ) , n = 1, 2, . 



Cn+l = (|A 



O n ,£ 2 O r 



-1/2 



(35) 



where Ao = 0. Since J x is hermitian, and CJ X is antiher- 
mitian, all O n can be chosen to be hermitian (antihermi- 
tian) for even (odd ) n, since A„ are purely imaginary. 
Thus it is easy to prove, for both inner products defined 
in (|23[) and (|2~4")) . that C has no diagonal matrix elements 

(o ni Cdn^J = 0. Also, it is straightforward to prove by 



induction that \O n ) is an orthonormal set: 

(d n ,d n >^ = S nn >, n,n' = 0, 1, . . .oo. (36) 
In this basis, C is given by the Liouvillian matrix, 



L, 



(O n \£\O n , 



Ai 

_ | Al A 2 

A* A 3 I 

AS ... 
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In both T = 0, oo limits, the current fluctuations func- 
tion (|30p is given by the root expectation value of the 
resolvent: 



G(z) = {j x ,{z-L)- x J x ) 



(38) 



where fia = (J x , J x ) is the corresponding zeroth moment. 
In these limits, all moments of Eq. (|33|) are root expec- 
tation values: 



Mfe — (Jx, £ Jx] 



(39) 



Using (|37p . an explicit, and useful relation between re- 
currents and moments is obtained 



Hk[A] = ^ {L k ) 
|Ar| 2 , 
|Ai| 4 + 



oo 



/'2 
/'4 
A*6 



+ |Ai| 2 |A 2 | 2 , 
|A 1 | 6 + 2|A 1 | 4 |A 2 | 2 
+|A 1 | 2 |A 2 | 2 |A 3 | 2 , 



|A 1 | 2 |A 2 



(40) 



D. Continued fraction representation 



Inverting z — L in Eq. (j3"5|) using elementar y alg ebra 
yields the continued fraction representation [38|, [39| 



G(z) = ^o- 



|Aj| 2 

. ia 2 i : 



(41) 



/io, are infinite. This form amounts to replacing the de- 
nominator's self energy by a purely imaginary constant, 



F(z) = 



1 



z - Ei (2) 
|A!| 2 

AS 



(44) 



(ii) Constant recurrents, A„ = A, n = 1, 2, . . ., yield 
semicircle imaginary part, 

F(z) 



z — 



A 2 



F» = v/1 - (o;/2A)2. (45) 

(iii) Linearly increasing recurrents, |A„| 2 = nfl 2 /2 define 
the continued fractions 



F(z) = 



1 



(46) 



which is relevant to HCB at high temperatures, as argued 
in Section IIV Al In Appendix [B] it is shown that ImF is 
a Gaussian, 



F» = ,%exp l-- 2 



(47) 



Incidentally, the real part of F(lu) is the Dawson func- 
tion [Zol 



(48) 



At zero temperature we obtain, 

G T ^{z) = 2(0|J 2 |0) , (42) 

Z ~T^ — 

where |0) is the ground state of H . Similarly, in the high 
temperature limit, the leading order G(z) is given by 

G T =oo{z) = ^Tr (Jl) | Ar|2 ■ (43) 

Z |Ag°| 2 
v s 



The infinite list of recurrents fully determines G"(u>). If 
only a finite set is computable, extrapolation of |A„| 2 to 
large n is unavoidable. Some intuition about the relation 
between recurrents and the imaginary part of the con- 
tinued fraction function at z —> to + it is gained by the 
following examples, for complex functions F(z), where 
F(uj + ie) = F'\<jj)-iF"{<jj). 

(i) A Lorcntzian, as given by the Drude form (JTJ) , is a 
somewhat 'pathological limit: All its even moments except 



E. Finite temperature corrections 

At finite temperatures, the current fluctuations function 
includes the effects of the thermal density matrix 

p{(3) = 2 N e^ H /Z(P). (49) 

(The prefactor of 2 N is introduced for later convenience) . 
One can write 

G(fi,z) = {{ P ,Jx},Q{z)Jx) 00 

n 

G n (z) = ^(z-L)-},. (50) 

Taking z — > u + it yields an orthogonal polynomial ex- 
pansion for the current fluctuations function, 

00 

G"(P,u>) = J2Cn{0)Pn(u>)G o (u;). (51) 

n=0 
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The polynomials P n arc orthogonal under the measure 
defined by Gq{uj), 

poo 

/ (kjP n (u))P m (u).G ((j) = £„ >m (52) 

J oo 

In Appendix we derive explicit expressions for P n as a 
function of Go and the preceeding recurrents |A m | 2 , m = 
1,2,. ..n. 

IV. KUBO FORMULA FOR HCB AT HALF 
FILLING 

In the previous section we introduced the moments and 
recurrents of the current fluctuations function G(z). 
Henceforth, we specialize to HCB at half filling. First 
we analyze the expected asymptotic behavior of the Li- 
ouvillian matrix elements. Second, we construct a varia- 
tional harmonic oscillator (VHO) basis in which the cur- 
rent fluctuations can be expanded. Third, we generate a 
high temperature expansion of G"(/3,w). 

A. Liouvillian of HCB and Gaussian asymptotics 

The Liouvillian hyper-operator C (|25[) . describes strong 
scattering in the following sense: When C acts on a hy- 
pcrstate composed of n spins, 

A \ 1 Q!l,a2,...CKn Q&l QOC2 QCtn 

A n — / y c ii,j?.....i„. J ii °i 2 ■■■°in ' 

the number of spins increases or decreases by precisely 
one, i.e. 

£A n — A n _)_i 

+ (53) 

When |j4 n +i| >> the primary effect of £ is to 

proliferate the number of spins. 

Let us compare the behavior of HCB with weakly inter- 
acting bosons. Consider a typical boson liquid Hamilto- 
nian 

H weak = a t i j ofl + i 5a t oa t a , (54) 

The action of the respective Liouvillian on a linear oper- 
ator yields 

[H weak , a f ] = a^H + go) act. (55) 

Thus, when g is "smaller" than Hq, the primary effect of 
£ is the first term, which propagates a' , rather than the 
second term which proliferates it. 

The root hyperstate Oo ^ 4, is bilinear in spins. By 
repetitive applications of £™Oo one obtains clusters of up 
to n spins. Consider a lattice in two dimensions or higher 
[ill ], with coordination number z > 2. For a typical A n , 
there are nz e g available bonds to attach an extra spin 
to an existing cluster, where z e g < z. Therefore, the 
number of distinct terms in the resulting A n+ i is roughly 



a factor of nz c g more than the number of terms in A n . 
We thus expect, for such a lattice, most of the weight of 
the hyperstate O n to consist of n spin operators. 
Since O n , O n _i are normalized, we can crudely estimate 
the asymptotic n dependence of the recurrents 

|A„| 2 « (6 ni C6 n ^ 2 ~ {4J) 2 z eff n, n » 1. (56) 

(The factor 4J stems from the commutation relations 
[S a , Sp] = ie a fs-ySj and the definition of the Hamiltonian 
§5§ ). As we have seen in (pf?|). if the asymptotic relation 
(|56|) were precise, the continued fraction expansion for 
G" would lead to a perfect Gaussian of width 4^2z c sJ. 
However, as demonstrated in Fig. 0] the linearity of |A„| 2 
is not precise: When C acts on O n it generates some small 
fraction of n — 1 spin operators, which are not included 
in the term A n O„_i of Eq. ([55]) . This spoils the exact 
relation between the index n, and the number of spin 
operators. 

Nevertheless, the coefficients C n {j3) in Eq. (|51[) are ex- 
pected to converge rapidly at high temperature. Expand- 
ing C n (f3) in a high temperature series, 

C n (/3) = ]TC^, (57) 

i 

it can be verified that the contribution of any m-spin 
operator A m to order j3 l depends on traces such as 

Tr [{H\ J x }A m ] = 0, m>i + 2. (58) 

H) 

Therefore, Ck measures the relative weight of i + 2 spins 
operators or less, in O n . Since, as argued previously, 
these weights become smaller for large n, we can expect 
for fixed i that 

lim -» 0. (59) 

ri>>i 

Thus at a finite order (3 l , only a finite number of Cn^ 's 
are of substantial magnitude. 

Phrased differently: In an idealized situation in which 
O n would contain only products of (n + 2)-spins, the C n 's 
would decay with (3 as 

C n = iTr (e-P H {6 n , Oo}) ~ ^CW + o(/3"+ 2 ). (60) 

This discussion raises an obvious question: Is there a 
lattice where Eq. 160\) becomes precise ?. We have prelim- 
inary expectations [42|, that this would be the case, at 
least for low orders in n, for the Bcthc lattice in the limit 
of large coordination number. The rapid convergence of 
the harmonic oscillator basis discussed below, indicates 
that the square lattice at high temperatures is not "far" 
from the infinite dimensional limit. 



B. Variational harmonic oscillator expansion 

A high temperature expansion for the conductivity can 
be generated by choosing a convenient basis to expand 
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G"(ta). The linear increase of the recurrents, suggesl 
earlier in Eq. ([55]) implies a Gaussian decay of G"( 
at high frequencies. Therefore it is natural to choosi 
variational harmonic oscillator (VHO) basis such thai 

G"{fi,w) = J2 D n ^)H n {u)^(u), (( 

n=0 



where H n (cj) = H n (u> / 'fl v ) j \J1 n n\ , and 



1 / CJ 

exp 



2 



(< 



The function i/jq(u>) is the VHO ground state probabil 
density, with a variational frequency scale £l v , and H n { 
are the corresponding (normalized) Hermite polynomh 
which constitute an orthogonal set under the measure 

The moments of H n (uj)ipo(aj)' 2 are 

IT 

for k > n. By construction, all vanish for n > k. 
Using the high temperature expansion of a finite number 
of moments 



vho 




FIG. 2: Infinite temperature current fluctuations function 
G'^ui). Solid lines (blue online) vho depicts the results of 
the variational harmonic oscillator expansion, and (Green on- 
line) extr depicts the result of the reccurrents extrpolation, 
shown in Fig. [3] Circles (red online) ed depict the exact 
diagonalization result computed on a 4 x 4 lattice. Inset: 
Convergence of the variational harmonic oscillator expansion. 
| AG" | are the distances between consecutive approximants, 
and n m ax is the number of computed moments. 



Mfcl 



(P) = Y,KP t , fc = 0,2,...7w, (64) 



and similarly, 



(65) 



We insert (foTj) into (|33[) and obtain a finite set of linear 



i(0 



equations for each Dn 



(66) 



n=0 



Solving for (|66p yields the desired coefficients for Eq. 



V. DYNAMICAL CONDUCTIVITY: HIGH 
TEMPERATURE 

The dynamical conductivity, to leading order in j3, is 



(67) 



We have calculated the infinite temperature current fluc- 
tuation Gy =oc , using three distinct methods. The results 
of the three approaches show satisfactory agreement, as 



k 







4 


2 


64 


4 


4096 


6 


544768 


8 


1.20906 x 10 8 


10 


3.96113 x 10 10 


12 


1.75571 x 10 13 



TABLE I: Low order moments at T = oo. 



depicted in Figf2] First we computed the infinite tem- 
perature moments, listed in Table |T1 



^Tt{J x C k J x 



0,2, 



12. 



(68) 



These traces were computed numerically on a finite 16 
site square lattice with periodic boundary conditions. At 
k > 8 the calculation introduces finite size errors, due 
to loops of operators which circulate around the system. 
We eliminated most of the contributions of these loops by 
averaging over Aharonov-Bohm fluxes through two holes 
of the torus § . 

The VHO calculation of G" follows Section QVB] Invert- 
ing Eq. (|66p . for the n max lowest moments, we obtain 
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FIG. 3: The recurrents of G^(w) at infinite temperature, e 
their extrapolation to large index n. 



FIG. 4: Dynamical conductivity in the high temperature re- 
sistive phase. The function can be fit by Eq. (|73[) 



G" moi . We obtain the variational frequency 

Q v = 7.48 J, (69) 

which minimizes the distance between the two highest 
order approximants, 



\ag"\'(q v ) 



J —oo 



G'i 



f °° duj 

J —oo 



G' 



(70) 



We plot G" for n max — 12 as the solid (blue online) curve 
in Fig. [2] In the inset of Fig. [21 we plot the convergence 
as a function of n max . The rapid decay of AG" indicates 
convergence of the VHO expansion at T = oo. 
The second calculation extrapolates the recurrents, as 
shown in Fig. [3J We compute |A£°| 2 , n = 0,1,... 6, 
from the known moments, using (|40p . The higher order 
recurrents exhibit an approximate linear increase with 
n — 1 which is consistent with a Gaussian decay of G" 
at high frequencies as given by the continued fraction 
example of Eq. ([17)) . We continue the linear slope (see 
dashed line in Fig. [3]) by the extrapolation, 



|A~ >6 | 2 



z e// (4J) 2 (n-l), z e// = 3.59, (71) 



where z e f y is the "effective coordination number" . The 
approximate function G" extr is obtained using the con- 
tinued fraction representation (|43|) . and plotted as the 
solid (red online) curve in Fig. [21 We notice that the two 
approaches yield very similar curves. 
Last, we compute the infinite temperature current fluc- 
tuations function by exact diagonalization (ED) in the 
Lehmann representation, 



GM - % E IH + E n E„ 

n.m 



(72) 




f-sum rule, 
high T expansion 



kinetic energy 



FIG. 5: The temperature dependent kinetic energy (— K x ) 
(solid, red online), calculated by Refs. @, @, and the sum 
rule for the high temperature expansion of the dynamical con- 
ductivity of Fig. [4j (solid, blue online). The vertical dashed 
line marks the superconducting phase below Tbkt- 



where \n) , E n are the eigenstates and eigenspectrum of H 
on a 4 x 4 lattice with periodic boundary conditions. We 
expect the finite size effects to be small at infinite temper- 
atures, where the correlation length is much shorter than 
the lattice size. Indeed, the agreement between the three 
approaches shown in Fig. [21 supports this expectation. 
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A. Finite temperature corrections and f-sum rule 



resistivity as 



Order (3 2 corrections to G" are obtained by a high tem- 
perature expansion of the moments //j^ as defined in 
_____ ______ 

Eq. (|64[) . Inverting Eq. (f6"6"j) . the coefficients Z?n were 

computed up to order n max = 12, and inserted into Eq. 
(|65]) . The resulting temperature dependent dynamical 
conductivity is plotted in Fig. [4] for several temperatures. 
A crude analytical approximation is given by 

q 2 tanh(w/(2T)) 



a(7» 



0.91- 



4.8J. 



(73) 



The truncation errors in the VHO expansion, and the 
high temperature series, are monitored by comparing the 
integrated conductivity to the kinetic energy, using the 
f-sum rule equation 



dcu tanh(/3w/2) 



G"(p,u) = {-q 2 K x ) p . (74) 



The values of {K x )p were computed by high temperature 
series expansion up to 11th order in Ref. [4__j, and at low 
temperatures by quantum Monte-Carlo simulations (__], 
e]. In Fig E] we plot the f-sum rule, Eq. ([74]), with G" 
calculated to order (3 2 , and the temperature dependent 
kinetic energy. The corrections grow larger than 10% at 
T < 2Tbkt- 

We note, however, that numerical satisfaction of Eq. ([74"]) 
does not, in general, ensure the accuracy of the DC con- 
ductivity c(0). In fact, there is an expectation of a cusp 
at zero frequency arising from high order non-linear cou- 
pling between the current and other long lived diffusive 
modes [____§. The magnitude of this feature is expected 
to be weak because of particle hole symmetry. From our 
calculations we can conclude that it is below the resolu- 
tion provided by the exact diagonalization and the 12th 
order moment expansion 45 1. 



R HN = R Q - 



T 



(76) 



where D v is vortex diffusion constant, and n v is the den- 
sity of free vortices. Estimating that n v = £(T)~ 2 using 
(fT2")h the critical resistivity is expected to be suppressed 
toward Tbkt as 



Tbkt \ 
T — Tbkt J 



1/2N 



R HN (T) = i? Q ^I-U- 2 cxp \-2B 

' (77) 

where no is a microscopic density of order a 2 . 

The precise diffusion constant for the HCB model is not 

known, and can be expressed as 



D v ng = nJ/H, 



(78) 



where k is an undetermined dimensionless constant. We 
can estimate the numerical value of k by requiring that 
R HN (T) matches Eq. (75]) at T > 2T BKT where the f- 
sum rule is satisfied to a higher accuracy. We find that 
matching occurs in the range, 



0.38 < k < 0.51. 



(79) 



VI. DC RESISTIVITY: HIGH TEMPERATURE 



We can interpret the vortex diffusion constant D v = 1 2 /t, 
as arising from a scattering time of order t = h/ J (the 
inter-site vortex hopping time, according to calculations 
of Ref. [H) and a short mean free path I fa a. Eqs. 
(|T8"]) . (fT9"]) are consistent with the value D v w h/m which 
was posited [46| for helium films. 



VII. DYNAMICAL CONDUCTIVITY: ZERO 
TEMPERATURE 



At T = the conductivity is given by 

.oM^ 2 ^^)+ q r , (aj) - 



(80) 



Based on the calculations of cf(T,uj) of the previous sec- 
tion, the high temperature expansion of the DC resistiv- 
ity, shown in Fig. [I] is given to order 1/T as 

R{T) = 0.23 R Q - (1 - 2.9( J/T) 2 + 0(( J/T) 4 ), ) (75) 

where Rq = h/q 2 is the boson quantum of resistance. 
In the regime of T > 2Tbkt, the negative corrections of 
order T -3 are relatively small, as the sum rule shows in 
Fig.[5] 

As Tbkt is approached from above the resistivity drops 
rapidly. The critical regime was described by Halperin 
and Nelson [H, [H, Hi)] (HN) who considered the contri- 
bution of unbound vortices to the charge transport coef- 
ficients. Using vortex-charge duality, and Einstein's re- 
lation for vortex conductivity, HN derived the critical 



We calculate the current fluctuations G T=0 in two stages. 
First, we appeal to the relativistic Gross Pitaevskii field 
theory [48[ , and obtain the low frequency gap and thresh- 
old form of the function. The mass parameter m, and the 
high frequency Gausssian scale f2o are variational fitting 
parameters. Second, we compute the lowest 12 moments 
by exact diagonalization of a 20 site cluster. These de- 
termine the lowest recurrent parameters |A°| 2 which are 
fit to the variational form. 



A. Field theoretical calculation 

The relativistic Gross-Pitaevskii (RGP) model describes 
the long wavelength theory of quantum rotators in the 
absence of a linear time derivative term 1481 . We use it 
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FIG. 6: Low order current fluctuations function in the rel- 
ativistic Gross Pitaevskii field theory, as calculated by Eq. 
(|84|) . Dashed line describes gapped magnitude fluctuations. 
Solid line describes gapless phase fluctuations. 



FIG. 7: The recurrents of G't= (lj) a t zero temperature. 
Dashed line emphasize the different behavior of even and odd 
recurrents, which is related to the gap structure in Fig. [8] 



to describe the HCB action at half filling, which exhibits 

particle-hole symmetry. *(r,r) denotes the fluctuating where c * = V/V A2 is the s P eed of sound - We compute 
condensate order parameter, governed by the imaginary the integral by performing the Matsubara sum, and take 
time action q ^ and (3 — > oo , which yields 

Shop = j<?xj dr ||*|* + ^|W| 2 -^ (|*| 2 - A 2 ) 2 , G^H - £ (^~) |^ e (M - "0- ( 85 ) 

(81) 

where A is the ground state order parameter, and m is The dynamical conductivity therefore exhibits finite fre- 
the effective (Higgs) mass. Expanding © to second <l uenc y WG1 S ht above the mass § a P m - We cannot rulc 
order about * = A(l + r^e**, the harmonic fluctuations out that hl S her order interactions may produce subgap 
arc described by decoupled phase and magnitude modes, spectral weight, since magnitude excitations can decay 

into phason pairs [27] . 

( 2 ) i ( 2 f ('2 2 \2\ In addition, the RGP theory Eq. (|8Tj) only describes long 

S rgp = 2 / d x J dr ^ J wave length fluctuations. Therefore, lattice scale and 

f f / p \ nigh energy cut-off effects are important for the high frc- 

+| / d 2 x / drA 2 [j] 2 + ^-|V?y| 2 + mifj . quency tails of G" . These will be computed directly from 

the zero temperature Kubo formula. 

(82) 



We introduce a vector potential by shifting V — ► V+iqA. 
The current operator is obtained from J x = 8S/5A X , 
therefore its paramagnetic part, expanded to the same 
order as Eq. (|82|) . is given by 

J x (x) ^q Ps {\7 ( f> + 2i 1 \7 ( f ) ). (83) 

The lowest order current fluctuations function is given by 
the bubble diagram [ZiJ [5(| depicted in Fig. [HI 




K; 

(cs 2 (tv n + v n ) 2 + (q + k) 2 + Cs 2 m 2 ) (cs 2 v 2 + k 2 ) ' 

(84) 



B. Variational Fit of Recurrents 

The moments of G'^ =0 , Eq. |33|) . are equal to the ground 
state expectation values, 

Mfc = (0|{J x ,£ fe J a }|0). (86) 

We compute a set of Hk, k < 12 by exact diagonalizations 
on a 5 x 4 lattice. Using Eq. ([40)) we determine the set 
of recurrents |A„| 2 , n = 1, 2, ... 6, which are depicted in 
Fig. \7\ Finite size effects are not expected to be large 
for small n, where A„ depend mostly on short range cor- 
relations. We notice a striking difference between even 
and odd recurrents, which seem to follow two linearly in- 
creasing slopes. This even-odd effect of the recurrents, is 
an indicator (not a proof) of a gap-like structure [Hlj . 
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0.08 



0.06 



02 




FIG. 8: Zero temperature dynamical conductivity of HCB. 
The Higgs mass gap m is depicted. Solid line (blue online) 
is determined by variationally fitting the recurrents to those 
in Fig. [7] Dashed line (red online) is computed by exact 
diagonalization of a 20 site cluster, with 5-function broadening 
of 0.5 J. 



Motivated by the field theoretical calculation, we use the 
trial function, 



fin 

\JFrp_ 



G 



RGP 



(w) exp(-w 2 /fig) 



(87) 



where m - the mass gap, and f^o - the high freqeuncy fall 
off, are variational parameters. We use (|40|) to determine 
the trial recurrents. The variational parameters are then 
determined by a least squares fit between the exact and 
trial recurrents. 
Thus we obtain, 



A test of these calculations is provided by the zero tem- 
perature f-sum rule is, 



duj G't =0 (uj) 



n 2 . 
1 J0+ 



Ps 



(90) 



The left hand side for the varational and exact diago- 
naliozation results yields 



du> G't =0 (io) 



o+ 



0.0148 J Variational 
0.0164 J ED 



(91) 



The sum rules are comparable in magnitude to values 
obtained by QMC [f| for the zero temperature kinetic 
energy and superfluid stiffness, 



-K* 



Ps 
Ps 



1.09765(4) J, 
1.078(1) J, 
0.019(1) J. 



(92) 



We notice the very small spectral weight (2%) of the high 
frequency peak at zero temperature, relative to the con- 
densate weight. This weight is due to the quantum fluc- 
tuations of the ground state, and to the non conservation 
of the current operator in two dimensions. This weight 
can be ascribed to the magnitude oscillations mode. 
Note: An important question is whether the mass gap 
survives corrections to Eq. (|84| . Within our variational 
approach, we tried to answer this question by allowing 
sub-gap spectral weight parameterized by a power-law 
tail such as 



Grr 



T=0 



oc 



(M/m)° 
1 + (M/to) q 



exp(-w 2 /^). (93) 



The least squares fit of the recurrents has found that 
a tends to increase indefinitely. This is consistent, (al- 
though not being a proof) with having a true gap at 



ft -► 8.6J, 
m -> 10.9J. 



(88) 



The mass gap appears to be similar to the high fre- 
quency Gaussian fall-off of the conductivity at high tem- 
peratures, as given by the linear slope of the recurrents, 
Eq. CLE]). 

The resulting dynamical conductivity is depicted in Fig[5] 
The figure also includes, for comparison, the results of 
exact diagonalization on a 16 site cluster, as given, for 
uj > by 



G'M = nJ2 I <0| J x \m)\ 2 5{u; + E - E m ). 



(89) 



The oscillations in G" ed are artifacts of finite size gaps in 
the spectrum. We see that the two curves agree on the 
position of the central peak and the total spectral weight. 



VIII. DISCUSSION 



A. Bad Metallicity 



HCB at half filling exhibits "bad metal" characteristics, 
as demonstrated in Eq. ((2]) and in Fig. [TJ 
In contrast to conventional metals and bosonic gases at 
high temperatures [10, Hi], [22j , the resistivity of HCB 
rises approximately linearly, without a sign of saturation 
at R w h/q 2 . Such behavior signals the breakdown of 
Boltzmann equation, since the mean free path becomes 
shorter than interparticle distance [l9j ]. 
A related quantity is the width of the low frequency 
conductivity peak, which in metals is called the "Drudc 
peak" . Here, Eq. ^ shows that the low freqeuncy tem- 
perature dependent peak in cr(u>) is goyerened mostly by 
the fluctuations-dissipation factor [52|]. If one fits the 
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width by Eq. U) in the regime T B kt < T « Vl 
obtains, 



= 2T. 



We note however, that for our hamiltonian ([5]), Tq 
not much smaller than f2. Therefore, separating thi 
frequency scales in a(T, lo) is difficult. However, the 
of the two scales can made larger by additional int 
tions. 

The HCB model relates the asymptotic resistivity 
to the zero temperature supcrfluid stiffness: 



dR c 



dT 



0.245 



Rq 



In a three dimensional system of weakly coupled k 
the transition temperature T c is shifted from Tbkt 
factor given by Eq. (JT3J) 0, [S3| . Above T c the dens 
free vortices rises rapidly, which causes the rise in 1 
as shown in Fig.[T] We can operationally define a non 
critical normal state resistivity R c by 



R c 



dR° 
~~dT 



(96) 



Thus, a HCB version of "Homes law" is obtained: 

p s {0) = 0.245^T C . (97) 

tic 
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FIG. 9: Comparison of "Homes law" of hard core bosons 
(HCB) and cuprate supeerconductors. HCB calculation is 
depicted by the dashed line (blue online) . Cuprate data, com- 
piled from Refs. [24|, HU, describes the two dimensional critical 
conductivity a c , as measured by optical conductivity exprap- 
olated to zero frequency, at the transition temperature T c . 
p 3 values are compiled from measured plasma frequencies by 
Eq. (|100|l . Rq = 6453f2 is the boson quantum of resistance. 
Different data points of the same symbol correspond to dif- 
ferent doping concentrations of the same compound. 



B. Cuprate Conductivity 

A large linear in T resistivity [23| and optical relaxation 
rate [53l . [H| are widely observed in clean samples, espe- 
cially near optimal doping. It has been shown [55[ that 
the linear resistivity is not consistent with a proximity to 
a quantum critical point. 

It is plausible therefore, that the "bad metal" character- 
istics of the normal phase of cuprates may be described 
by lattice bosons in their resistive state. Support to 
this viewpoint is give n by Uemura's empirical scaling law 
T c cx Pg b (T = 0) pjll, and the observation of a superfluid 
density jump in ultrathin underdoped cuprate films [56| . 
These are consistent with the behavior of a bosonic su- 
pcrfluid, captured by an effective XY model. In under- 
doped cuprates, additional evidence exists that the hole 
pair bosons survive above T c , up to the pseudogap tem- 
perature scale HH T* > T c . Thus there have been 
several theoretical approaches [13, \E^, [62| to the 

superconducting properties of cuprates based on lattice 
bosons of charge q = 2e. 

We note however, that above the pseudogap tempera- 
ture and frequency scale T* , hole pairs completely dis- 
integrate. Thus, above T* , the HCB model cannot be 
used to obtain the temperature and frequency dependence 
correctly. 



Homes et al. [24(, have pointed out, that the superfluid 
stiffness is generally proportional, with a seemingly uni- 
versal constant, to the product of "critical conductivity" 
a^ D , times T c . The "critical conductivity" a c was ex- 
perimentally defined by extrapolating cr(u>, T c ) to zero 
frequency. This empirically universal scaling was inter- 
preted as a sign of " dirty superconductors" , where T c is 
determined by the disorder driven scattering rate [o3j . 
Here we promote a different viewpoint, which maintains 
that Homes law can also be obtained in a disorder free 
model. The necessary ingredients are strong scattering 
effects, as described by hard core bosons above Tbkt- 
To relate the experimental data of cuprates to the HCB 
model, we first translate the 3D critical conductivity to 
a two dimensional critical conductivity: 

a c = a c a 3 c D , (98) 

where a c w 1.5nm is the interplanc distance. 

The zero temperature boson superfluid stiffness, can be 

deduced from the in-plane London penetration depth, 

where c is the speed of light, or the optical plasma fre- 
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qucncy aj ps : 

h 2 

P g = — 2 acUJ p- ( 10 °) 
lone z F 

In Fig. [9] we plot the data reported in Ref. 0, b Y 
translating it into the two dimensional quantities. We 
find, quite remarkably, that the slope of Eq. (|97[) lies 
quite close to the experimental data. We emphasize that 
the proportionality constant in Eq. (|97|) is not expected 
to be universal. A priori, one might expect it to vary with 
additional interactions and the filling. If the agreement 
with cuprates' data is not fortuitous, it would imply that 
this constant might not be sensitive to moderate varia- 
tions of H. 

At very low temperatures, we find that there is a signa- 
ture of the quantum fluctuations of the order parameter 
in the optical conductivity at high frequencies. These are 
magnitude fluctuations characterized by a "Higgs mass" 
gap at frequency 

m « 10p s . (101) 

It is interesting to mention that such a massive mag- 
nitude mode, given by the relativistic Gross-Pitaevskii 
action, appears in the superfluid phase of the strongly 
interacting Bose Hubbard model at integer filling. It has 
been associated with the "oscillating superfluidity" ex- 
periments of cold atoms in an optical lattice [25|, |2g, [27j • 
While a clear signature of such a peak in cuprates has not 
been identified, we speculate that it might be related to 
the ubiquitous mid infrared peak which has been detected 
in several compounds at low temperatures [28[ . We cau- 
tion however, that at these high frequencies, additional 
fermionic excitations become increasingly important. 
In summary, we conclude that some of the "normal state" 
phenomenology of cuprates, and perhaps other uncon- 
ventional superconductors, may be described by lattice 
bosons. However, we emphasize that Eq. (JSJ) oversim- 
plifies these systems, by omitting potentially important 
ingredients: fermionic excitations, long range Coulomb 
interactions, intcrlayer coupling, disorder and inhomo- 
geneitics, and of course the HCB conductivity away from 
half filling. Clearly, these effects need to be accounted 
for in future work. 
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APPENDIX A: ORTHOGONAL POLYNOMIALS 
AND RECURRENTS 

Iterative use of the matrix inversion formula for G n (z) 

G n (z)=» (z-L)-) ) , (Al) 
with L defined as in (|37[) yields, for z — > lu + ie 

g;;m = p n ( w )G (w), (A2) 

where P n (uj) is a polynomial that depends on the |A m | 2 's 
for m < ri, and is given by 

n 

P n (w) = JJ|A fc |det (A3) 
fc=i 

In the above, [z — L) 1 is the upper left (n— 1) x (n— 1) 
sub-matrix of (z — L). 

The determinants det (z — L) n _ 1 obey the recursion re- 
lation 

dot (z - L) n+1 = z det {z - L) n + |A„| 2 det (z - L)„_j . 

(A4) 

Equation (|A4[) is a recursion relation for orthogonal poly- 
nomials. The polynomials P n {uj) defined in Eq. (|A3[) are 
therefore orthogonal polynomials under the scalar prod- 
uct defined by 

/oo 
n in • 
(A5) 
-oo 

The complexity of computation of C n ((3), P n {u}) depends 
on obtaining all the low A rn ,m < n. Therefore, the 
expansion (|5ip would be useful if it could be truncated 
at finite n, provided that the coefficients C n ((3) decay 
rapidly enough with n. 

APPENDIX B: LINEAR RECURRENTS AND 
THE GAUSSIAN SPECTRAL DENSITY 

Let us consider a sequence of recurrents |A„| 2 given by 

|A„| 2 = i7ir! 2 n=l,2,... (Bl) 

The continued fraction representation can be solved using 
the following map. Consider the dimensionlcss position 
operator x represented by raising and lowering operators 
of the one dimensional harmonic oscillator 

.x=^(a t + a). (B2) 



Since we have (n + l\x\n) = y/(n + l)/2, the function 
G n {z) of Eq. (f5TJ|) is equivalent to 

G n (z) = ^{n\—^—\Q). (B3) 
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Using the x representation for the ground state of the where N n = 1/ y/2 n nl. 
harmonic oscillator we have 

f°° 1 e~ x2 



cuj 2 /Q 2 

W fi IH / e* dt 



lu + ie — Qx ^/tt 

1 + - / 

(B4) 

Likewise, for higher values of n we have 

T r ( M t f°° j N n H n (x) e-* 2 
ImG n (w)//i = -Im / dx — — =- 

J-oo (J + le-ilX y/TT 

(B5) 
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